BUHEP-06-03 



Dark Matter in the Simplest Little Higgs Model 
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We investigate a variation of the Simplest Little Higgs (SLH) model that is invariant under a 
discrete Z2 symmetry. Our motivation for imposing this symmetry is the two dark matter candidates 
it implies, a scalar -q and a heavy neutrino. We examine the viability of these two candidates using 
the standard relic abundance calculation. Direct detection prospects are also discussed. 



I. Introduction 

The need for Dark Matter (DM) , confirmed by current 
astronomical observations, is a sufficient reason to look 
beyond the Standard Model (SM). A reliable way to in- 
clude DM into a model of physics beyond the Standard 
Model (SM) is to introduce a discrete symmetry under 
which some new particles are odd, while all Standard 
Model (SM) particles are even. The lightest odd particle 
is stable and is a candidate for DM. In supersymmetry 
(SUSY) this discrete symmetry is i?-parity, in Universal 
Extra Dimension (UED) models it is KK parity, and for 
Little Higgs models this discrete symmetry is called T 
parity 0,1,0. 

In this paper we will study the DM candidates of a T 
parity invariant version of the Simplest Little Higgs 0, Q 
model. Unlike other T parity invariant Little Higgs mod- 
els 'g', , it is not possible to make all new particles odd 
under T parity. The gauge boson Z' remains even. Since 
tree level Z' exchange is allowed, the Z' needs to be heavy 
to avoid conflict with precision electroweak experiments. 
This means we must take the scale / to be relatively high. 
The price we pay for this is fine tuning. Despite this fine 
tuning, the lightest T parity odd field, which is either a 
scalar 77 or a heavy neutrino remains as an interesting 
DM candidate. The scalar 77 is similar to the electroweak 
singlet limit of the LH scalar DM found in |^]. It anni- 
hilates predominantly into higgs boson and tt pairs, but 
the cross section is small unless 77 is quite heavy. The 
heavy neutrino is a new addition to Little Higgs model 
DM. It is also distinct from SUSY neutralino and UED 
heavy neutrino DM since it is nearly inert under the SM 
weak gauge group. The SLH heavy neutrino does couple 
to the heavy gauge bosons and to the higgs. 

We investigate the viability of both of these DM can- 
didates using the standard relic abundance calculation. 
As an additional check, we also explore the constraints 
direct detection experiments place on SLH DM. In an ef- 
fort to make this analysis more general, we do not limit 
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ourselves to a particular value of the symmetry break- 
ing scale /. Hopefully this will allow our analysis to be 
applicable to other models similar to SLH. 

We find that 77 can be an acceptable DM candidate 
only when coannihilating with a slightly more massive 
fermion. The 77 mass that yields the proper abundance 
will depend on the scale / and the type of fermion it 
coannihilates with, but typically ~ 1.0 TeV. We find 
77 — nucleon elastic cross section to be well below the 
current experimental sensitivity for the 77 mass range of 
interest. The SLH heavy neutrino is an acceptable DM 
particle over a wider region of parameter space. How- 
ever, direct detection constraints are also much stricter. 
Current direct detection experiments rule out scenarios 
with / ^ 3.0 TeV, and are predicted to be sensitive to 
/ ~ 4.5 TeV by the time they finish running. To be con- 
sistent with the above bounds, heavy neutrino DM must 
have TOjv ^ 600 GeV for / = 3.0 TeV and tojv ^ 1.0 TeV 
for / = 4.5 TeV. 

This paper is organized as follows: In section II we in- 
troduce the SLH model and show how it must be modified 
to make it Z2 invariant. We describe some of the exist- 
ing constraints on this model from precision electroweak 
observables in section III. After a brief overview of the 
standard relic abundance calculation in section IV, we 
examine the DM candidates 77 (section V) and the heavy 
neutrino (section VI) in more detail. In section VII we 
include the effects of coannihilation, which we find to be 
substantial. Section VIII is devoted to the experimen- 
tal constraints on SLH model DM from direct detection 
experiments. 



II. Fields and Representations 

In the Simplest Little Higgs model, the SM weak gauge 
group is extended to SU{3)u]'S)U{l)x- AH SM electroweak 
doublets are therefore accompanied by a new partner to 
form SU triplets. We follow the second set of repre- 
sentations given in SLH since it has no SU{3)^ anomaly. 
In this set of representations, the first and second gener- 
ation quarks are promoted to a 3 of SU{3)u], while the 
third generation quarks and the leptons are a 3. All fields 
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and their representations under {SU{3)c, SU{3)w)u(^i)^ 
are listed below. 



,c3 ,,c3 



*Q3 = (3,3)j 

(3,1)- 



*Q1 

d 



*Q2 = (3,3)o ^'Lz = (l,3)_i 

u-\u-^ = {3,l)^ < = (l,l)o 
$i,$2 = (1,3)^ 



Under this assignment, the partner of the third gener- 
ation SM quark doublet is an up type quark T, while 
the first and second generation partners are down type 
quarks D, S. The partner of the lepton doublet for each 
generation i is a neutrino Ni. 

This extended weak gauge group is broken down to 
SU{2)w^U{1)y when two different scalar SU{3) triplet 
fields $1, 3>2 acquire vacuum expectation values. 
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Of the 10 pseudo-Goldstone boson (PGB) degrees of free- 
dom in $1 and $2, 5 are eaten by the gauge bosons of 
(5t/(3),„ (g) U{l)x)/{SU{2)w «) U{1)y). The remaining 5 
PGBs decompose under SU{2)w (Xi U{1)y as a complex 
doublet h, our higgs, and a real scalar rj. 

The PGBs feci a potential generated radiatively by 
gauge and Yukawa interactions. However, no single cou- 
pling breaks enough of the global symmetry protecting 
the higgs to give it a mass. Terms that include more 
than one coupling can give the higgs a mass, but these 
terms are higher order in <i>i,$2 and the higgs mass 
terms they generate depend only logarithmically on the 
cutoff scale A. This 'collective breaking' is a necessary 
ingredient in all Little Higgs models 0] . Along with 
gauge and Yukawa interactions, this model also contains 
a tree level scalar potential V = /i^$|$2 + h.c. This 
potential gives a mass to 77 and is needed to achieve the 
correct higgs vev with a higgs that sufhciently heavy to 
avoid experimental bounds. 

T parity operator 57 

We now want to impose a Z2 symmetry under which 
SM particles are even and new particles are odd. The 
new particles predicted in this model are the fermions 
T,S,D,N,, the scalar 77, and the S'C/(3)/S'C/(2) gauge 
bosons 

We begin by looking at the SU{'i)w fermion triplets. 
We define the action of T parity on the fermion triplets 
to be 



Q,L, 



(2) 



where Vt — Diagonal{^l, —1, 1). As desired, the first two 
components of ^'i, a SM electroweak doublet, are even, 
while the third component is odd. By defining ^ with 



a minus sign, we have ft € SU{3). This simplifies our T 
parity implementation greatly. Since our Lagrangian is 
manifestly SU{3) invariant, all factors of Cl are automat- 
ically taken care of. All that remains to be determined 
is a ± sign for each field. We can determine the sign by 
whether the field is a SM field or a new field. For the 
higgs triplets, the transformation is 
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(3) 



which requires /2 = /i = f/V2- Under this transfor- 
mation, the higgs doublets are T parity even, while the 
remaining PGB, 77, is T parity odd. The vevs of the <I> 
are also T parity even. 

Having defined the action of T parity on SU{3) triplets, 
we immediately know how the octet of gauge fields A^j, 
transform. 
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(4) 



This transformation leaves the diagonal gauge bosons 
and the gauge bosons in the upper left 2 by 2 block 
even, while all other gauge bosons are odd. This is 
immediately a problem since the Z', a new particle, 
remains even. For now we press on with our discrete 
transformation, but we will discuss the consequences of 
an even Z' in more detail in the next section. 

Enforce T parity on the Yukawa sector: 

To determine how the rest of the fields transform, we 
need to enforce T parity invariance on the Yukawa sector. 
To be consistent with the T parity transformation of the 
higgs, this means that both $1 and $2 must couple with 
equal strength. 

For the first two generations of up quarks, the bottom 
quark, and the leptons, enforcing T parity is straight for- 
ward and requires no additional operators. The fields 
rf'^'^, u'^^, m"^^, and are all T parity even. Enforcing T 
parity on the down quark Yukawas is slightly more diffi- 
cult, for reasons that become clear by examining the top 
Yukawa in greater detail. 

The operator that leads to the top mass is 



t„,c3\ 



(5) 



The transformation properties of the u^^ are deter- 
mined by requiring that the linear combination of the 
u^i^ that gets a mass with the heavy top quark is odd. 
Expanding out the higgs fields to lowest order, we find 



-^^{uf - uf)Q, + K^{uf + uf)T + 



V2^ 



(6) 



The correct T parity transformation for the uf^ is there- 
fore 
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(7) 



Comparing © to the corresponding equation in SLH, 
we see that we have halved the number of parameters: 
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/i, /2, A„2 Au,/- The reduced number of param- 
eters leads us to a simple relation between the mass of 
the top quark and the mass of new heavy quark T, 
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mt 



V2f 



(8) 



For / = 2.0 TeV this ratio gives m-r = 2.0 TeV, almost 
twice its value at the 'golden point' of SLH. 

The problem with the first and second generation down 
type Yukawas in SLH is now clear. The operator 

Ad*Qi($idf + $2^2') (9) 

has exactly the same form as (0 , and would imply a rela- 
tionship similar to JHJ between the SM down quark mass 
and the heavy D quark mass. For nid ~ 10 MeV This 
relation predicts a heavy down quark with a mass less 
than a GeV! Clearly we need to use a different operator 
to give mass to the heavy down type quarks. 

To avoid this catastrophe, we introduce operators that 
give mass to the T parity even and T parity odd portions 
separately 



+ $2)^^ + ed*Qi($i - $2)dL , (10) 

where — {d'i^+d2^), d'^ — (d j^ — d^^) are respectively 
the odd and even combinations [33 of The new 

particles and SM quark masses are no longer related for 

The only remaining term is the mass term for the heavy 
neutrinos Nij-n'i, both of which should be T parity odd. 
The operator 



A„n^($I + $t)^^ 



(11) 



is sufficient. 

The couplings (|10lll|) violate the little higgs mecha- 
nism. By this we mean that they generate the operator 
^1^2 + h.c. at one loop, and therefore lead to quadrati- 
cally divergent radiative higgs mass terms 
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(12) 

The larger the mass of the heavy quark or neutrino, the 
larger the contribution to the higgs mass. Since we are 
not aiming to remove fine tuning in this model, we will 
not restrict particle masses on account of (|12|) . However 
it does encourage us to look for light ^1.0 TeV DM. 

Now that we have made the SLH model consistent with 
T parity we can identify DM candidates. The DM par- 
ticle will be the lightest T parity odd particle (LTOP). 
The neutral T parity odd fields available are the neu- 
trinos (iVi, n^), the scalar 77, and the heavy gauge bosons 
Q . The mass of the heavy neutrino can easily be a few 
hundred GeV for a small value of A„, and we expect the 
mass of r/ (~ /i) to be approximately the weak scale. The 
Wp Q obtain a mass mw' — ^ when SU[?))w ® U{l)x 
breaks to SU{2)w ® U{^)y- The W mass and all of its 
interactions are fixed for a given scale /. Because of this 
lack of ficxibility and because the W is usually heavier 
than the ry and heavy neutrino, we will not consider the 
W' as a DM candidate here. 



III. Parameters and Tunings 

As we are unable to make the Z' odd under T parity, 
we need it to be very heavy in order to avoid conflict with 
precision electroweak observables (EWPO). The only 
way to get a heavy Z' in the SLH model is to increase the 
symmetry breaking scale /. The lower bound on / was 
originally estimated to be / ^ 2.0 TeV jjj, but later anal- 
ysis point to a stricter bound of / > 4.5 TeV p^lTlIT^. 
Rather than abide by one of these bounds, we will inves- 
tigate SLH DM for a variable scale / > 2.0 TeV. This 
allows our analysis to be more general, and applicable to 
other similar models. Whenever results for a fixed / are 
desired, we will use the value / = 4.0 TeV. 

For a given /, the masses of the heavy gauge bosons 
and heavy top are automatically determined. 
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The radiative contributions of these heavy particles to 
the higgs mass can then be determined by the usual tech- 
nique as shown in [li- As a consequence of a high value 
of /, the masses in (fT^ can be quite large, causing the 
theory to be somewhat finely tuned. However, since it is 
not our goal in this paper to remove fine tuning, we will 
accept it at its current level. 

For relatively light DM (^ 1 TeV), wc must also check 
to make sure that any effective four fermion operators 
are appropriately suppressed. Four fermions operators 
come about in this scheme through box diagrams with 
SM external legs and T parity odd particles on the inter- 
nal loop. The constraints on new contributions to four 
electron operators (eeee) are the most stringent, followed 
by constraints on four quark operators (qqqq). We have 
estimated the contributions of heavy particle loops to 
(eeee) , (qqqq) in the model we have presented, and we 
find that these contributions are less than 1/(5—10 TeV)^ 
in the DM mass region of interest. The (qqqq) contribu- 
tions from D — T] boxes are suppressed by factors of [v/ f) 
at each vertex (basically we have decoupled the heavy 
quarks from their SM partners by introducing the extra 
operators). Charge conservation forbids a eNrj vertex, so 
only W' — N box diagrams contribute to (eeee). These 
contributions are safely suppressed by the larger mass of 
the W. 



IV. Relic DM Abundance Calculation 

We now give a brief overview of the relic DM calcula- 
tion before examining the heavy neutrino and 77 further. 
We follow the standard procedure outlined in ll4L Il5l| . 

The number density n of a cold dark matter particle x 
obeys the Boltzmann equation 
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(14) 
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Here H is the expansion rate of the universe, Ueq is the 
equiUbrium number density, and uVrei is the annihilation 
cross section of the particle times the relative velocity. 
The ( ) around the annihilation cross section indicates 
that we take a thermal average. We will use the standard 
value for the equilibrium density of a cold, nonrelativistic 
particle of mass m with g degrees of freedom: 



/mT\3/2 



= 5 



V 27r / 



exp{-'m/T). 



(15) 



Since we expect the cold DM particles to be slow, we 
first take the non-relativistic limit of the cross section. 
To obtain the nonrelativistic cross section, we substitute 
s — 4M^ + M^w^ into the cross section and keep only 
first order terms in . Thermal averaging then gives 



(JVrel ~NR 0- " 
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X 



(16) 



where x = 

From the thermally averaged cross section we can de- 
termine the contribution of particle x to the total energy 
density of the universe. The density of particle x divided 
by the critical density corresponding to a flat universe, 
denoted as fi^/i^, is determined to be 



1.04 X 10^ 



Xp 



M, 



pi 



3b 



(17) 



to consider particles with approximately the same mass 
as the LTOP we must drop this assumption and examine 
the evolution of the total number density of T parity odd 
particles — n^^ + n-^^ + • • • . The evolution equation 
for the total number density takes the familiar form 



— + 3Hn 
at 



(19) 



Here CTe// is the sum of all cross sections cry — <j{xiXj ~* 
XX') weighted by the degrees of freedom gi in Xi and the 



mass difference = {m^^ 
hilating particles, we have 



J ) / . For N coanni- 



N 



a^jf = ^ ayif^(l + A,)3/2(l + A^.)3/2^-.(A.+A,)^ 

^i 9eff 



where 



9eff 



-^.g,(l + A,) 



3/2, 



(20) 



(21) 



Note that the mass appearing in a; = above is the 
LTOP mass. The effects of particles much heavier than 
the LTOP are exponentially suppressed, so only particles 
very close in mass to the LTOP are relevant in the sum 

®- 

The solution to H19|l is the same as (|17I18|I but with 
modified a and b. 



Here Xi;^ is the critical temperature below which the ex- 
pansion term alone determines the evolution of the num- 
ber density. It is referred to as the freezeout temperature 
and it can be expressed analytically as 



Xp 



log( 



0.047 g Mpi (a + |^) 



(18) 



The parameter g* is the total (spin, color etc.) num- 
ber of relativistic degrees of freedom at temperature xp- 
Typically, xp k. 20 

The results (|17ll8f) are derived by solving the 
Boltzmann equation assuming constant entropy per 
comoving volume in the limits x <^ xp and x 3> a^F, 
then matching the two solutions together. A more 
cornplete discussion of this procedure can be found 
in [13L [T5II . We will consider DM to be cosmologically 
acceptable if falls within the 2 a limits from 

WMAP, 0.094 < ^^h^ < 0.129 



where 



Xp 



1.04 X 10^ 



Xp 



pi 



dx , Ih 



(22) 



2x p 



dx 



The averaged coefficients ag/ / , 6e/ / are defined by sub- 
stituting aij,bij for ffij in H20|l . 

We now examine the relic abundance calculations for 77 
and N in more detail to locate the regions in parameter 
space where they are suitable DM candidates. For a given 
value of /, we will vary /x, An,,and (or equivalently 
TO^, TO^r, mu). For simplicity we will first consider the 
scenario where the LTOP is much lighter than the other 
new particles. Once we understand the abundances in 
the rriLTOP mj^p limit, we will explore the effects of 
coannihilation. 



Coannihilation 

If there is a particle {X2) .just slightly heavier than the 
LTOP (xi)j then coannihilation processes X1X2 ~^ XX', 
etc. become important in determining the LTOP num- 
ber density 0|. To study the effects of coannihilation 
we must modify our relic abundance formalism slightly. 
The formalism of H17I I18|) assumes that all T parity odd 
particles have decayed into the LTOP by xp. In order 



V. Scalar r? as the LTOP 

The first dark matter candidate that we consider is rj, 
the remaining PGB from Sf7(3)^ (g) U{l)x ^ SU{2)w 
U{1)y breaking. It receives a mass ^ /i^ from the 
scalar potential. 

To determine the viability of i] dark matter we need 
to calculate the annihilation cross section. As rj is an 
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EW singlet, it does not couple directly to the or 
gauge bosons. The r\r\Z' couplings is also zero because 
77 is a real scalar field. The dominant annihilation mode 
for m,, ^ / is into a pair of higgs bosons 7777 — *■ ft-o^o- 
This annihilation proceeds through the rj^h'^h term in the 
scalar potential, and also through the i-channel exchange 
of a heavy Wq gauge boson. 

The other important annihilation mode for i] is 7777 — > 
tt. As the rjTt'' coupling is large, we might expect 
7777 — *■ tt through T exchange to be the dominant annihi- 
lation mode. However, this process is helicity suppressed, 
meaning that this process has a = (see Eq. I16f) . The b 
term for 7777 — > ti can be large, but its effect on the over- 
all cross section is suppressed by an additional factor of 
a;F~O(20). 

There is also a contribution to 7777 — s- tt from higgs ex- 
change. The interaction responsible for this annihilation 

comes from the scalar potential and has strength ~ 
For heavy 7/ this becomes 0(1). Annihilation modes of 
77 to light fermions or SM gauge bosons are all negligi- 
ble, suppressed by small fermion masses or by additional 
powers of {v/f). 

In Figure 1. we plot the regions where fljjh^ is allowed 
by cosmology as a function of the mass of rj and F. We 
assume a higgs mass (3^ of 140 GeV. For a given /, we 
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FIG. 1: Regions in the rrir, — f plane where Qrjh^ is within 2 a 
of the WMAP limits. In this scenario we assume all other T 
parity odd particles are much heavier than 77. 

see that fth^ does not fall into the allowed range until 
the rj mass is 0{f). Because of the higgs exchange con- 
tribution, 77 annihilation into tt continues to grow as 777^ 
is increased. When niri ~ /, the cross section is elevated 
to the point that we get the right relic abundance. Once 
rriri is much larger than / the cross section is too big and 
77 DM is insufficient. Although there is a 771^ — / region 
where fi^/i^ is within the WMAP experimental limits, we 
cannot accept it as a viable DM candidate. The rj mass 
which gives the right relic abundance is always heavier 



than the mass of the W , and thus the 77 in this scenario 
would never be the LTOP. As we will see in section VII., 
the acceptable 777,, — / region changes significantly when 
we allow coannihilation with nearly degenerate fermion. 



VI. Heavy Neutrino(s) as the LTOP 

The heavy neutrino is a slightly more complicated 
LTOP since it is a dirac fermion. To calculate its anni- 
hilation cross section we must calculate the annihilation 
cross sections of each of it's Weyl components {N,n'') 
with themselves and with each other. 

at^t = cr{NN XX') + airfn" XX') + 
a{NN XX') + airfff -> XX') 
+ aiNn" XX') + aiNn" XX') + /i.e., 

(23) 

where X, X' are some SM particles. Though this looks 
tedious, many of the sub-cross sections are either zero 
or very small. The only significant terms are a{NN 
XX') and af^Nn" XX'). 

Unlike many heavy neutrino DM candidates [Tsl IT^ , 
the SLH heavy neutrino interacts very weakly with SM 
gauge bosons. It does not couple to the W^, and 
its coupling to the Z° is suppressed by a factor of 
gp-(l — tan^ dw)- AH annihilations proceed through the 
exchange of a heavy gauge boson, a heavy fermion, or 
the higgs. 

The NN annihilation cross section contains the mode 
NN i/D, £~£~^ through i-channel W exchange. For 
light 77iAr this is the dominant mode as the coupling of 
the N to the W is large. For heavier rriN, NN an- 
nihilation to fermions via s-channel Z' quickly becomes 
the most important process. The reason for this domi- 
nance is simply the presence of a pole at 771 = Mz' /2 
in the s-channel propagator. In vicinity of itin = Mz> /2 
the cross section is so large that its nonrelativistic limit 
must be treated specially. The usual Taylor expansion 
about V — (Eq. (llBfl ) is very inaccurate and can result 
in negative cross sections |3Jj. To avoid this inaccuracy, 
we follow one of the prescriptions suggested in and 
set u = in all s channel Z' propagators. 

For rriN 3> Mz'/'2,, Nn'^ annihilation into heavy SM 
quarks via higgs exchange may also become relevant as 
the hNn'^ coupling is proportional to ^^ijj^. It is in- 
teresting to note that, contrary to many DM candi- 
dates ia,ii|ll|23,|2l|, the annihilation into W'^W is 
not very significant for either DM candidate considered 
here. 

We have calculated the relevant terms in H23|l and the 
resulting Vl^qh^ is plotted in Figure 2. for both fixed / 
and variable /. 

In the 777 AT — / plane, we see two strips where the il^h? 
is consistent with experiment. The 777 values of the 
both strips increases with /, yet they remain roughly 
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n,(GeV) 

FIG. 2: In the top plot: We show the regions where Q.Nh^ is 
within the 2 a WMAP limits for a range of / and mjv. In 
the bottom plot: Q.Nh'^ for a heavy neutrino LTOP versus 
the mass of the heavy neutrino in GeV for fixed / = 4.0 TeV. 
In both plots we have assumed that all other T parity odd 
particles have mass ^ mjv. 



400 GeV apart. Between the two strips, heavy neutrino 
DM is unacceptable because there is not enough of it, 
while outside the strips there is too much DM. We can 
understand the origin of these two regions by looking 
at the relic abundance for fixed / as a function of the 
heavy neutrino mass, shown in the rightmost plot above. 



Right at the Z' pole, the annihilation cross section blows 
up causing the relic abundance to plummet. However, 
as mjv moves away from the resonant value Mz' /2 by 
±200 GeV, VtNh? crosses the allowed WMAP region. 
The acceptable m^v regions always border the Z' pole. 
Thus as / is increased, raising Mz' , the allowed neutrino 
DM mass also increases. 



VII. Including Coannihilation 

As we saw in section IV, particles slightly heavier than 
the LTOP can have a significant effect on the relic DM 
abundance. The nature of this effect depends on the 
self-annihilation cross section and the number of degrees 
of freedom of the heavy particle, as well the annihilation 
cross section of the heavier particle with the LTOP. 
Coannihilation of the LTOP with a particle that has the 
same (or fewer) degrees of freedom and participates in 
the same interactions usually results in a smaller cross 
section, while coannihilation with a strongly interacting 
particle with many degrees of freedom can increase the 
cross section by an order of magnitude or more. In the 
Simplest Little Higgs model there can be coannihilation 
among more than one generation of heavy neutrino, 
between the heavy neutrinos and the heavy quarks, and 
also between rj and a heavy fermionjs^. 

Heavy neutrino LTOP including coannihilation 

We first consider coannihilation among three degener- 
ate heavy neutrino flavors. Two neutrinos of the same 
flavor can annihilate with each other as described pre- 
viously, and two neutrinos of different flavor can annihi- 
late through the i-channel exchange of a W' . The flavor 
changing processes NiNj are less efficient 

than the self-annihilation processes, and as a result the 
average cross section at a given mass m^r for three neu- 
trino flavors is smaller than the cross section for a single 
neutrino flavor. This results in a larger iljv/i^, as can be 
seen in Figure 3. 

In order to achieve the correct relic abundance with 
three neutrinos, we need a larger self annihilation cross 
section to compensate for the extra degrees of freedom. 
This requires mAr to be even closer to Mz'/2, and thus 
the two acceptable DM regions are closer together. 

Coannihilation of a heavy neutrino with a heavy quark 
has a very different effect. Since they are colored, heavy 
quarks can annihilate to SM quarks through gluon ex- 
change, and thus a{DD XX') and cr{d%djj XX') 
include O(a^) terms. The quark self-annihilation cross 
sections are consequently orders of magnitude larger than 
a{NN XX'), especially for small nfiMTTno- As the 
heavy quarks couple to the Z', their cross section is also 
enhanced near the pole at = Mz'/'2. 

The large heavy quark self annihilation cross sections 
dominate the average cross section cr^ff since they are 
weighted by the large number of degrees of freedom in 
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n,(GeV) 

FIG. 3; Top: Acceptable QmH'^ for a single flavor of heavy 
neutrino (black) and including coannihilation among three 
flavors of heavy neutrino (gray, green online) versus mjv, /. 
All neutrinos are degenerate. Right: fijv/i^ for one (solid 
line) and three (dash-dot line) degenerate neutrino flavors for 
/ = 4.0 TeV. Dashed lines are WMAP hmits. 



the colored heavy quarks. The mixed annihilation chan- 
nel a{ND — > XX') through W' exchange is enhanced 
by color factors relative to the mixed neutrino annihila- 
tion cross section, but it is still much smaller than the 
self annihilation cross sections a{NN XX'), a{DD 
XX'). However, because aeff is so strongly controlled by 



a{DD -> XX')_and a{d%d% XX'), the inefhciency 
of the mixed ND processes has a negligible effect on the 
average cross section except for when mD,'mN are very 
close to mz' /2. 

In Figure 4. we plot the acceptable regions of heavy 
neutrino DM including coannihilation with a heavy 
quark as a function of tojv, /■ If the heavy quark is just 
slightly heavier than the heavy neutrino, we see that 
there can be an additional neutrino mass region where 
heavy neutrino DM is cosmologically allowed. The range 
of •mul'Tfijs! where this second region occurs is very 
small. For mD/'rnM > 1.05, this additional allowed mass 
region either does not exist or exists where is exper- 
imentally forbidden. Even if 1.0 < mu/m^ < 1.05 , the 
additional region is only acceptable when / > 2.5 TeV. 
For lower /, heavy quark annihilation through gluons is 
so dominant that (Jeff is increased to the point where 
there is insufficient DM. 

Coannihilation of rj with a heavy fermion: 

Following the same procedure as above, we can calcu- 
late Qjjh^ including coannihilation with a heavy fermion 
(either heavy neutrino or heavy quark). The coannihi- 
lation channels rjN XX' and r/Z? ^ XX' are medi- 
ated by heavy neutral W exchange. The heavy fermions 
have larger annihilation cross sections and more degrees 
of freedom than rj, so their self-annihilation controls cTe//. 

To show the effect of a degenerate heavy fermion on 
Q.rfh? , we plot Qj^h"^ including coannihilation with a heavy 
neutrino in Figure 5. 

The most important consequence of coannihilation in 
this scenario is that the s channel Z' exchange terms 
originally in the neutrino self annihilation cross section 
are folded into to the effective cross section for ry. As a 
result there is now a large dip in Jl^/i^ near M2//2. For 
toat within 10% of this dip is substantial enough that 
we achieve acceptable UnIt-^ for much lighter m^. Com- 
paring figures JQ) and (O , we can see that coannihilation 
decreases the allowed rurj by almost a factor of 10. If to^v 
is even closer to the dip becomes larger and we get 
the correct relic abundance at lighter m^. 

We now summarize the results of our relic abundance 
calculations before investigating possible experimental 
signatures. In order for flrjh^ to be consistent with the 
region allowed by WMAP for an 77 that is actually the 
LTOP, we must have coannihilation with a nearly de- 
generate fermion. Coannihilating with a heavy neutrino 
with mAf/m^ = 1.1 for / = 4.0 TeV, we achieved the 
correct relic abundance for ~ 1.0 TeV. If the LTOP 
is a heavy neutrino, f^jv/i^ falls within the experimental 
limits for a much wider range of mjv and /. Taking / 
fixed and varying m^r, there are two regions of acceptable 
V.jyh'^ - one on each side of the pole in the cross section at 
mN — Mz' /2. For / = 4.0 TeV the two allowed neutrino 
mass regions are rnjv ~ 900 GeV and rriM ^ 1300 GeV. 
If the other heavy neutrinos are nearly degenerate with 
the LTOP, the allowed mass regions at a given / are 
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FIG. 4: Left: Regions of cosmologically allowed Q,Nh^ in- 
cluding coannihilation with a heavy down type quark. Three 
different quark masses are shown; mD/mN = 1-1 (gray, green 
online), 1.05 (light gray, yellow online), degenerate (solid 
black). See figure ^ for Q.Nh^ without coannihilation. As 
before the allowed range for ^Nh^ is 0.094 < f7iv/i^ < 0.129. 
Right: Q,Nh^ versus mjv including coannihilation with a de- 
generate heavy quark for / = 4.0 TeV. The dashed lines 
indicate the WMAP limits. 



closer to Mz'/2. Including coannihilation with a heavy 
quark, we can achieve adequate values for flffh'^ for mjv 
as low as 200 GeV. To achieve a tolerable ftf^h^ at such 
a low m^v, mu/mN must fall in a very specific range. 




M^(GeV) 

FIG. 5: Regions of acceptable fi^/i^ including coannihilation 
with a heavy down neutrino, versus m^, /. Three different 
N — i] mass differences are shown; Am — mN — mn = 10% rriN 
(black), 5% mjv (gray, green online), degenerate (light gray, 
yellow online). There is an additional region not shown here 
because it requires m,, > M^i/'. The allowed range for fijv/i^ 
is 0.094 < STiiv/i^ < 0.129. 



VIII. Possible Signatures of SLH Dark Matter 

Now that we have identified a few candidate scenarios 
for DM in the SLH model with T parity, we look to see 
what signatures each predicts. 

By looking for anomalous atomic recoil in systems of 
very cold atoms, an upper limit can be placed on the DM- 
nucleon elastic scattering cross section. Currently the 
most stringent limits come using ''"^Ge in the CDMS 
experiment located in the Soudan mine in Minnesota. 
CDMS places separate limits on the spin independent 
(SI) and spin dependent (SD) interactions of DM with 
nuclei. Spin independent interactions are generically pro- 
portional to the reduced mass of the DM - nuclei system, 
while spin dependent interactions are proportional to the 
spin of the nuclei. For larger nuclei {A > 40) and non- 
Majorana DM, spin independent interactions are usually 
orders of r nag nitude larger than the spin dependent in- 
teractions [Tg[l2^. 

Based off of the total DM-nucleon elastic cross section, 
we can also make an estimate of the rate for indirect 
detection. DM particles are collected over time in a 
massive stellar object like the sun and they annihilate 
with each other. The energetic neutrino remnants from 
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these DM annihilations reach the earth and can be seen 
in high energy neutrino detectors. This rate does depend 
somewhat on the dominant annihilation mode of the 
DM and other model dependent parameters [Tgl . 

Detection of heavy neutrino DM 

Since the heavy neutrino is a dirac fermion its coupling 
to the Z' has {V—A) structure. The vector portion, when 
combined with the vector part of a qqZ' interaction, con- 
tributes to the spin independent heavy neutrino-nucleon 
scattering cross section. One might expect this contribu- 
tion to be small given the weakness of the coupling and 
the large mass of the Z'. A calculation of the A^— nucleon 
effective cross section shows that this is not the case. The 
Z' couples to both components of a weak SU{2) dou- 
blet with equal strength. The effective proton/neutron- 
Z' couplings are then approximately three times the size 
of the quark-Z' coupling. The coherent interaction of a 
Z' with the protons and neutrons in a large nucleus can 
thus be sizable Following [H [Hill Hi we cal- 
culate the neutrino-nucleus scattering cross section (per 
nucleon) in ^'^Ge to be 

/ 2 Te V N 4 

fTw-„„c = 3.7 X 10-42 (^^j cm2. (24) 

Normalized according to the convention in [H , the 
cross section is nearly independent of the heavy neutrino 
mass. The only mjv dependence is in the nuclear form 
factors, which we will ignore here. The cross section 
(|24|l therefore applies to both neutrino mass regions for 
a given /. For / £ 3.0 TeV, the aN-n y . r fr om is 
excluded by the current CDMS data [IJllll. Between 
/ = 3.0 TeV and / = 4.5 TeV, (Jn-uuc is beneath the 
current CDMS limits, but within the predicted final 
sensitivities of the current run, CDMSII. Heavy neutrino 
DM in this region would certainly be visible in later 
phases of CDMS, if not earlier. In Fig. 6 we show how the 
A^-nucleon scattering cross section per nucleon compares 
with the current and predicted CDMS sensitivities [l^ 
as a function of /. 

There is also a SD neutrino-nucleus cross section from 
the axial-vector portion of Z' exchange. We estimate it 
to be several orders of magnitude smaller than the SI 
cross section and will therefore ignore its effects. For a 
total {SI + SD) nuclear cross section of O{10~^ ph) (per 
nucleon), we expect 0(1) event/year per km^ |2J| in a 
neutrino telescope. 

Detection of rj DM 

By the same procedure, we can calculate the effective 
77— nucleon cross section. Since 77 is a scalar, there is no 
spin dependent interaction. The full 77-nucleon cross sec- 
tion is just the spin-independent 7y-nucleon cross section. 

An 77 scatters off of a light quark (m, d, s) predomi- 
nantly through t channel higgs exchange. Other rjq rjq 




10' id" 

WIMP Mass [GeV] 



FIG. 6: The solid line indicates the current CDMS limit on the 
spin independent DM-nucleon cross section per nucleon as a 
function of the mass of the DM particle. In the shaded region 
we shown (JN-nuc wherever fijv/i^ is cosmologically allowed, 
starting with / = 2.0 TeV. We assume mN <^ m_D,m^. As / 
is increased we move down and to the right within the shaded 
region. The dotted line is the projected final sensitivity for 
CDMS-II. The dashed line shows the expected sensitivity of 
the next phase of CDMS, SuperCDMS A [g^. 



processes, either from s channel heavy quark exchange 
or from higher order Yukawa terms like rj^{q'^q), are 
small. The higgs exchange interactions are proportional 

-p- and to the quark masses, thus the strange quark con- 
tribution is the largest. To get the full SI 77- nucleon 
cross section, we must also include 77-gluon scattering 
a{'qg ~* rjg). This occurs through higgs exchange with 
a quark loop which emits two gluons. Loops containing 
heavy SM quarks (c, 6,i) or new heavy quarks {S,D,T) 
can make substantial contributions to C7(j]g rjg) since 
both as and the heavy quark Yukawa couplings are 
large |0,|23,'2^. Using a standard approximation for the 
heavy quark (both SM and T-odd) loops jSlJ we calculate 
the 77— nucleon cross section in the limit mT,mD ^ jti,, 
to be 

4^/ rn„ \2/2.0TeV\4 „ , , 

ar,-nuc 10 { — I cm^. (25) 

^ V500 GeYJ \ f J ^ ' 

This cross section is well below the current and projected 
CDMS limits for the entire mrj — f range of interest. 
Given cr,,_„„c ^ 10~^ pb and the lack of any spin 
dependent interactions with the nucleus, the potential 
for indirect detection of 77 in a neutrino telescope is also 
very dim [2^. 



IX. Conclusion 

One consequence of enforcing a 22 symmetry onto a 
Little Higgs model is that the lightest odd particle is a 
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potential DM candidate. We have investigated this here 
in the context of the Simplest Little Higgs model. We are 
unable to make all new particles odd, and therefore we 
are forced to work within a somewhat finely tuned model. 
Accepting this fine tuning, we examined the two DM can- 
didates in this model; a heavy neutrino {N, n'^) and a 
scalar 77. Through the standard relic abundance calcula- 
tions we have found the circumstances under which these 
DM candidates are allowed by cosmology. We have also 
checked to make sure these circumstances are consistent 
with current direct detection limits. 

The first DM candidate we investigated was the scalar 
•q. We found that 77 cannot be the LTOP unless there 
is a nearly degenerate (within 10%) fermion. Coannihi- 
lating with a nearly degenerate heavy fermion, there is 
parameter space where Vt^ih? is within the WMAP limits 
and is in the TeV range. For / — 4.0 TeV, this oc- 
curs for rriri ~ 1.0 TeV, although the actual number will 
depend on the type of fermion and the degree of degen- 
eracy. We estimate the 77 — nucleon cross section to be 
^ 0(10"^) pb in the entire region of interest. This cross 
section is well below the predicted sensitivity bounds of 
both the current CDMS run and the first SuperCDMS 
phase. ?7 DM in this scenario would be difficult to find. 

The other DM candidate we considered was a heavy 
neutrino. When the mass of the neutrino is close Mz' /2, 
the annihilation cross section is enhanced and lowers 
ripfh^ into an acceptable range. Typically this happens 
for \mN - Mz'/2\ ~ 200 GeV. Provided that the neu- 
trino mass is in the right range, flpfh^ stays in the al- 
lowed region for a large range of the other T-odd particle 



masses. However, heavy neutrino DM is ruled out by 
direct detection unless / > 3.0 TeV. Based on the pro- 
jected final CDMS II sensitivity (see fig.®), the bound 
on / could become as high as 4.5 TeV if no heavy neu- 
trino DM is seen by the end of the current CDMS run. 
These bounds are approximately independent, al- 
though CDMS constraints on more massive particles are 
somewhat weaker. The bound / ^ 4.5 is consistent with 
the strictest EWPO bound [13, Ell- Heavy neutrino DM 
for / = 4.5 TeV has mass ~ 1 TeV, although it can be 
much lighter if a heavy quark is nearly degenerate with 
the heavy neutrino. As neither the heavy neutrino nor 
r] has significant SD interactions with nuclei, we expect 
indirect detection signals of SLH DM to be very small. 

Some other distinct features of this model are the T- 
even Z' , a very heavy mr ''-^ f top quark partner, and 
the lack of heavy T parity even quarks In addition 

to the direct detection signals we discussed, this model 
and its DM candidates may yield interesting collider phe- 
nomenology which remains to be investigated. 



Acknowledgments 

First and foremost we would like to thank Martin 
Schmaltz for his help and guidance throughout this work. 
We are also grateful to Tuhin Roy for many helpful dis- 
cussions and for reading early drafts. This work was sup- 
ported by the U.S. Department of Energy under Grant 
No. DE-FG02-91ER40676. 



[1] H.-C. Cheng and I. Low, "TeV symmetry and the little 

hierarchy problem," JHEP 09 (2003) 051, 

[hep-ph/0308199 
[2] H.-C. Cheng and I. Low, "Little hierarchy, little 

Higgses, and a little symmetry," JHEP 08 (2004) 061, 

[hep-ph/0405243 
[3] I. Low, "T parity and the littlest Higgs," JHEP 10 

(2004) 067, hep-ph/0409025 
[4] M. Schmaltz, "The simplest little Higgs," JHEP 08 

(2004) 056, hep-ph/0407143 
[5] D. E. Kaplan and M. Schmaltz, "The little Higgs from a 

simple group," JHEP 10 (2003) 039, hep-ph/0302049 
[6] A. Birkedal-Hansen and J. G. Wacker, "Scalar dark 

matter from theory space," Phys. Rev. D69 (2004) 

065022, hep-ph/0306161 
[7] J. Hubisz and P. Meade, "Phenomenology of the littlest 

Higgs with T-parity," Phys. Rev. D71 (2005) 035016, 
hep-ph/0411264 
[8] M. Schmaltz and D. Tucker-Smith, "Little Higgs 

review," hep-ph/0502182 
[9] M. Perelstein, "Little Higgs models and their 

phenomenology," hep-ph/0512128 
[10] Z. Han and W. Skiba, "Little Higgs models and 

electro weak me asurements," Phys. Rev. D72 (2005) 

035005, 1 hep-ph/ 0506206 



[11] G. Marandella, C. Schappacher, and A. Strumia, 

"Little-Higgs corrections to precision data after LEP2," 
Phys. Rev. D72 (2005) 035014, hep-ph/0502096 

[12] J. A. Casas, J. R. Espinosa, and I. Hidalgo, 
"Implications for new physics from fine-tuning 
arguments. II: Little Higgs models," JHEP 03 (2005) 
038, hep-ph/0502066 

[13] M. Srednicki, R. Watkins, and K. A. Olive, 

"Calculations of relic densities in the early universe," 
Nucl. Phys. B310 (1988) 693. 

[14] K. Griest, "Cross sections, relic abundance, and 

detection rates for neutralino dark matter," Phys. Rev. 
D38 (1988) 2357. 

[15] G. Servant and T. M. P. Tait, "Is the lightest 

Kaluza-Klein particle a viable dark matter candidate?," 
Nucl. Phys. B650 (2003) 391-419, hep-ph/0206071 

[16] WMAP Collaboration, D. N. Spergel et. al, "First 
Year Wilkinson Microwave Anisotropy Probe (WMAP) 
Observations: Determination of Cosmological 
Parameters," Astrophys. J. Suppl. 148 (2003) 175, 
astro-ph/0302209 

[17] C. L. Bennett et. al, "First Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Preliminary 
Maps an d Basic Results," A strophys. J. Suppl. 148 
(2003) 1, 1 astro-ph/0302207| 



11 



[18] K. Griest and D. Seckel, "Three exceptions in the 

calculation of relic abundances," Phys. Rev. D43 (1991) 
3191-3203. 

[19] G. Jungman, M. Kamionkowski, and K. Griest, 

"Supersymmetric dark matter," Phys. Kept. 267 (1996) 
195-373, hep-ph/9506380 

[20] H.-C. Cheng, J. L. Feng, and K. T. Matchev, 

"Kaluza-Klein dark matter," Phys. Rev. Lett. 89 (2002) 
211301, hep-ph/0207125 

[21] R. Mahbubani and L. Senatore, "The m inimal mod el 
for dark matter and unification," hep-ph/0510064' 

[22] CDMS Collaboration, M. S. Armel-Funkhouser et. al, 
"Exclusion limits on the WIMP nucleon cross-section 
from the first run of the Cryogenic Dark Matter Search 
in the Soudan underground lab," Phys. Rev. D72 
(2005) 052009, astro-ph/0507190 

[23] M. Kamionkowski, K. Griest, G. Jungman, and 

B. Sadoulet, "Model independent comparison of direct 
versus indirect detection of supersymmetric dark 
matter," Phys. Rev. Lett. 74 (1995) 5174-5177, 
hep-ph/9412213 

[24] F. Halzen and D. Hooper, "Prospects for detecting dark 
matter with neutrino telescopes in light of recent results 
from direct detection experiments," hep-ph/0510048 

[25] M. Drees and M. M. Nojiri, "New contributions to 
coherent neutralino - nucleus scattering," Phys. Rev. 
D47 (1993) 4226-4232, hep-ph/9210272 

[26] M. Drees and M. Nojiri, "Neutralino - nucleon 

scattering revisited," Phys. Rev. D48 (1993) 3483-3501, 



hep-ph/9307208 
[27] G. Servant and T. M. P. Tait, "Elastic scattering and 

direct detection of Kaluza-Klein dark matter," New J. 

Phys. 4 (2002) 99, hep-ph/0209262 
[28] CDMS Collaboration, D. S. Akcrib et. al, "Limits on 

spin-independent WIMP nucleon interactions from the 

two-tower run of the Cryogenic Dark Matter Search," 

astro-ph/0509259 
[29] CDMS-II Collaboration, P. L. Brink et. al, "Beyond 

the CDMS-II dark matter search: SuperCDMS," 

ECONF C041213 (2004) 2529, astro-ph/0503583 
[30] All CDMS plots were generated using the following 

website: A Dark Matter Limit Plot Generator 

http://dmtools.berkeley.edu/limitplots, by Rick 

Gaitskell and Vuc Mandic. 
[31] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, 

"REMARKS ON HIGGS - BOSON INTERACTIONS 

WITH NUCLEONS," Phys. Lett. B78 (1978) 443. 
[32] The choice of which combination of (F^ to be odd/even 

here is arbitrary. We chose these particular 

combinations to be consistent with the u'f' . 
[33] Calculations with = 115 GeV and = 200 GeV 

yield almost identical results. 
[34] Strictly speaking, this inaccuracy is problematic only 

for TUN > r^/. We calculate Vz' ~ 70 GeV <^mN. 
[35] There can also be coannihilation of either DM candidate 

with the W' but we have neglected these effects here. 
[36] There is also a small contribution from Z exchange, 

which we include in all calculations. 



